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1 Introduction 


The cross-sections of QCD multiparticle production processes at high energies have many 
similarites with the multiparticle distributions of a classical gas, an analogy which was first 
noted by Feynman and Wilson |^. This gas is essentially one dimensional in rapidity space. 
In this paper we use the gas analogy to derive a partition function for the Lund string 
fragmentation model . We perform a virial expansion to the second order in the density of 
particles. Our partition function then yields an equation of state for a Van der Waal’s gas. 
Furthermore, it reduces to that of an ideal gas when the produced particles are massless. 

The partition function of the gas is related in a simple way to the multiplicity distribution 
of its constituent particles. This provides us with a method of investigating the partition 
function. We show that for a fixed rapidity “volume” our partition function corresponds to 
a multiplicity distribution which is very similar to a binomial distribution. 

For large rapidity intervals the major fluctuations in multiplicity stem from gluon radi¬ 
ation. We will, however, neglect gluon emission. In this paper we are only interested in 
comparing the Lund fragmentation model with the properties of a classical gas. 

We analyse the multiplicity distributions obtained from the JETSET Monte Carlo 
for several scenarios. Firstly, we investigate the string break-up vertices, then the primary 
particles and finally we include decays. We find that all cases are remarkably well described 
by distributions from the binomial family. In the derivation of our partition function we 
assume that the particles are ordered in rapidity. Since this is true for the vertices, we 
expect the distributions of vertices to be the optimal case. Indeed, these distributions are 
well described by our partition function. 

The transition from vertices to particles introduces some smearing in rapidity. This results 
in a wider multiplicity distribution, where the width is sensitive to the transverse mass of 
the produced particles. We obtain an ordinary binomial for the primary particles. However, 
the strong smearing from decays ensures that, for the final state particles, this distribution 
becomes a negative binomial distribution. 

We shall begin with a short presentation of the basic ideas of the Feynman-Wilson 
gas (FWG). This is followed by an introduction to the Lund model and its relationship 
to the FWG. We next turn to the multiplicity distributions for the vertices and lastly how 
they are modified for the final state particles. 


2 The Feynman—Wilson gas 

The original discussion of the FWG can be found in Q. Here we summarize the main features 
of the model. We consider a multiparticle production process where the two primary particles 
have four momenta pi and p 2 and large invariant s = {pi + P 2 )‘^- The n secondary particles 
have four momenta fci, k 2 , ■ ■ ■, kn, and each is on the mass shell. In the FWG model the 
three remaining degrees of freedom in each ki correspond to the “spatial” co-ordinates of a 
gas particle via 


X = kx 


1 


y = ky 


5 = 

where the transverse mass is dehned by 

\n[{kz + ko)/mA] 

(1) 

mj_ = 

Jw? + kl + k"^ . 

(2) 


Note that in this picture i corresponds to the rapidity of the relevant particle. We will 
assume here that each produced particle is of the same type (each has the same mass) but 
the extension to different species is straightforward. 

We can write the total cross section for the production process using these spatial variables. 
We hrst note that the invariant phase space d^k/ko becomes d^r. The energy momentum 
conserving delta functions are first written in terms of p = pi + p 2 — ki — ... — kn- 

S{Po)d^{p) =‘26{p+)6{p-)6‘^{p±) (3) 

with pj- = po ± Pz. This can be expressed in terms of r variables using the relationship 
ko± kz = 

The delta functions have the effect of introducing a hxed volume for the gas. The transverse 
momenta are limited and constrain the gas to a narrow tube of radius ~ 300 MeV. We shall 
instead focus on the f co-ordinate. We first introduce and W- via 

W± = {pi+p2)± ( 4 ) 

so that we can write 

<5(p±) = -^m_iiexp{±Zi)) . (5) 

In the following we use the Lorentz frame where W± = ^/s. The two delta distributions 
contain the requirement that the the “gas volume” should be of the order of Ins. To see this 
we may integrate out the rapidities of the first and the last particles to obtain 

dzidZnS{- ■ ■)(5(- ■ ■) ~ 1/s 

^ -Zn ^ In(ys) . (6) 

We may in this approximation choose a number sq in such a way that 

Az = zi- Zn = \n{s/so) ( 7 ) 


and assume that all the particles are kept inside this rapidity “volume”. If the spatial co¬ 
ordinates of the primary particles are Ri and R 2 respectively then the cross section can be 
written as 

For fixed Ri, R 2 then (7^ corresponds, in the FWG analogy, to the partition function of the 
gas and the functions are the n particle distribution functions for the gas. Our aim is to 
connect these ideas to particle production within QCD as represented by the Lund model. 


25{p+)d{p-)6'^{p±) ... ,rn,i?i,i?2) 


( 8 ) 
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The Lund model and the Feynman—Wilson gas 


3.1 The Lund model 

In this section we briefly review some features of the Lund model fragmentation scheme. 
We will mostly be concerned with the simple situation when the colour force held from an 
original quark-antiquark pair (produced by e’''e“ annihilation, for example) decays into a set 
of final state hadrons. 

In the Lund model, the colour force held is approximated by a massless relativistic string 
with a quark (q) and an antiquark (q) at the endpoints. The gluons are treated as internal 
excitations on the string held. This means that there is a constant force held, k ~ 1 GeV/fm, 
corresponding to a linearly rising potential, spanned between the original pair. After being 
produced the q and the q are moving apart and the energy in the held can be used to produce 
new qq-pairs. When a new pair is created the string is split into two pieces. 

The production rate of a pair with combined internal quantum numbers corresponding to 
the vacuum is, from quantum mechanical tunneling in a constant force held, given by 

= exp • (9) 


Here the quarks in the pair have transverse mass = \J mass /r and transverse 

momentum The final state mesons in the Lund model correspond to isolated string 

pieces containing a q from one breakup vertex and a q from the adjacent vertex together 
with the produced transverse momentum and the held energy in between. The break-up of 
the string is illustrated in Fig.(|^). 

One necessary requirement is that to obtain real positive (transverse) masses all the ver¬ 
tices must have spacelike difference vectors. Together with Lorentz invariance this means 
that all the vertices in the production process must be treated in the same way @]. Another 
consequence is that it is always the slowest mesons that are produced first in any Lorentz 
frame (corresponding to the fact that time-ordering is frame dependent). Furthermore each 
vertex has the property that it will divide the event into two causally disconnected jets, the 
mesons produced along the string held to the right and those produced to the left of the 
vertex. This can be seen in Fig.(p. 

A convenient ordering along the force held of the produced particles is rank ordering. 
Two particles have adjacent rank if they share a qq pair created at a vertex . The first rank 
meson contains the internal quantum numbers of the original q together with those of the q 
produced at the vertex closest to the endpoint q. Similarly the second rank meson contains 
the internal quantum numbers of the q from this “first” vertex and the q of the “second” 
etc. In this way rank ordering corresponds to an ordering along a light-cone. Alternatively 
it is also possible to rank order in the direction from the original q. 

The basic Lund model fragmentation process then stems from the following two assump¬ 
tions 
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Figure 1: The break-up in space-time of a Lund string into n hadrons. The fragmentation area is 
denoted by A. 



Figure 2: The production, in energy-momentum space, of a particle with transverse mass m_i_. 
The particle is produced between the vertices 1, with the squared proper-time rf = Ti/k^, and 2, 
with r| = T 2 /k^- The particle has fractional light-cone components z_|_ and z_. 

1. In the centre of phase space (i.e. far from the endpoints) the string decay process will 
reach a steady state. The probability to hud a vertex is, after many production steps 
along the light-cone, a hnite distribution in the proper time of the vertex. This is also 
the case when the total string held energy becomes very large. 

2. The decay process is the same whether it is ordered along the positive or along the 
negative light-cone. 

If we consider Fig.(^, this means that we assume that the probability to reach the space- 
time point 1 at (a;+i,a;_i), after many steps along the positive light-cone, and to produce a 
meson with transverse mass mj_ by one further step to the vertex 2 at (x+ 2 ,a:- 2 ), is equal 
to the probability to reach the point 2, after many steps along the negative light-cone, and 
by one further step to 1 produce the meson with m^. 

Changing variables to the squared Lorenz invariant proper time = x+x_ and the rapid¬ 
ity y = l/21n(x+/x_) the probability to reach the point 1 is HirDdrldyi. The probability 
to produce one further particle with mass m± and fractional light-cone component Z-\. is 
m±)dz+. A particle with fractional light-cone component has the positive light-cone 


4 





Figure 3: The decay, in energy-momentum space, of an n-particle cluster with invariant squared 
mass s. The fragmentation area of the cluster is Arest- T = is with respect to the proper-time 
r of the last vertex. 


energy-momentum component p_|_ = z+nx+i and has, in order to stay on the mass-shell, the 
negative component p_ = = Z-KX- 2 - This means that we obtain the equation: 

H{T‘^)dT^f{z+,m±)dz+ = H{T 2 )dT 2 f{z-,m±)dz- . (10) 


It is a nice and surprising feature of the assumptions above that there is a unique process 
that fulhlls Eq. (p!oD 

Hj = exp(-6r) with F = , 

exp{—bm\/z) . (11) 



The numbers Cj and Njk are normalisation constants and the particle is assumed to be 
produced in a step from a vertex with flavour j to a vertex with flavour k. If u/ denotes 
the number of qq-flavours, the process has Uf + 1 parameters. Although the parameter a 
is, in principle, flavour dependent, there has been no need to utilize this in the Lund model 
as implemented in the JETSET Monte Carlo; except for the hrst rank particle in a heavy 
quark jet [Q. The parameter b must be flavour independent. 

It is possible to construct the probability to produce a hnite energy cluster of rank- 
connected particles from Eq. ([TT|) . Such a cluster is shown in Fig.(^). This probability 
distribution is in a natural way subdivided into two parts, the probability to obtain the 
cluster and the probability that the cluster decays in a particular way. In the following we 
order the particles along the positive light-cone. If the cluster has a total light-cone fraction 
z and a fixed total squared cms energy s then the (non-normalised) probability to obtain 
such a cluster is 


dPext = 

Z 


z \ 1 — z 

exp(—6F(s, z)) with F(s, 2 ;) = s- 


( 12 ) 
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The cluster then starts at a vertex with flavour /q and ends with flavour /„. The T value 
is that of the last vertex, as shown in Fig.(||). Thus a cluster is produced in the same way 
as a single particle between the vertices with oq and a„. Similarly we And that the (non- 
normalised) probability for the cluster to decay into the particular channel with the particles 

{v]j is 


dPint = [JJ Njdpj5{p] -m]) di^Pj - Ptot) exp(-6Arest) 


(13) 


where A^est is the decay area of the cluster, as shown in Fig.(^). Equation (|^ is for simplicity 
written in the ordinary Lund model fashion with a single a-parameter (this parameter is not 
explicit in the formula) and we note the appearance of the phase space for the final state 
particles multiplied by the exponential area decay law. The quantity Ptot is the total energy 
momentum of the cluster so that = s. We may determine the finite energy version of 
the vertex distribution, P(F), from Eq. ([T^) by exchanging for F. This yields 


The function Hg in Eq.(P^ is exponentially decreasing in F so that the power dependence 
in the denominator only plays a role for small values of F and then it is hardly noticable 
for large values of s. In this way the assumption 1. above is fulfilled. That is to say when s 
becomes very large there is (after normalisation) a finite distribution in the proper-time size 
of the decay vertices. 


3.2 The connection between the Lund model and the FWG 


We will now exhibit the decay distribution of a cluster, as given by Eq. (0), in terms of the 
partition function which is studied in statistical physics. For simplicity we write the formulas 
for a single particle transverse mass m±_ and a single flavour and we let j denote the rank 
of a particle. The phase space factor can in analogy with the result in section ^ be written 
with the particle energy momentum vectors Pj = mj_(exp( 2 /j), exp(—j/^)) as 


dT = [^Ndpj5{p] - m]) 5(^pj - Ptot) 


n 


1 

t2 ^-1 


(5(^ mx exp(|/j) - P+)6{Y^ exp(-j/j) - P_ 


N- T-r 
— 

2 


3 ■ 


(15) 


We have in the last line integrated out the first and the last rapidities in the delta function 
and from now on we assume that the remaining particles are placed in rapidities between 
Ay/2 and —Ay/2 with Ay = ln(s/so) and Sq is some suitable scale. If all the particles are 
ordered in rapidity we may integrate out the phase space factor and obtain 


/ 




7V2(iVA|/)^-2 
s{n — 2)! 


(16) 
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Figure 4: The fragmentation area partitioned into two-particle regions reveals how the area can 
be expressed as in Eq.(lT^). 


We next consider the decay area of the cluster. Figure shows that it can be written in 
terms of the rapidities of the particles 

n n 


A = ml^Y^ exp(|/fc - yj) . 

j=i k=j 


(17) 


From this equation and Eq. (lT^) we note that the decay distribution in Eq.(|I^) has similarities 
with a partition function, and we therefore dehne a grand partition function Z as 


n n 

- 


Ndyj I (5(.. .)5(...) exp ( —hm\ EE exp(2/fc - yj) 




j=i k=j 


n \ / ^ n n 

Y[Ndyj] 6{...)6{...)exp 

j=l / \ j=l k=j 


(18) 


(The factor of s is required in order to have a dimensionless partition function.) In this way 
we see that the decay distribution in Eq. (|1^) may be interpreted as the partition function for 
a system of n particles with co-ordinates yj interacting with exponential two-body potentials 
in a one-dimensional volume equal to Ay. We note that whilst all the particles interact in 
this way (“long-range interactions”) the exponential decrease of the potentials ensures that 
the effective interaction is rather short ranged. 

If the particles are imagined as making up a gas in rapidity space and are interacting via 
two-body potentials, then the Hamiltonian is 


j j,k 


(19) 


The phase space volume element is with Hj denoting the quantities canonically 

conjugate to the co-ordinates yj. The kinetic energy factors T are integrated out in Eq.([T8|) 
and incorporated into the constants N. These constants then play the role of fugacities. 

We shall now attempt to obtain a simplifed expression for the partition function, We 
have already seen that in the strong ordering limit the phase space may be easily integrated 
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according to Eq.(^). Using this limit may seem drastic since two neighbours in rank may 
well have a different rapidity order. However, if many pairs are not well ordered then many 
of the exponential potentials in the in the partition function will be strongly increasing, i.e. 
the area suppression in the Lund model will make these contributions small. 


We can easily find an expression for the exponential in our partition function if we approx¬ 
imate the fragmentation area. Assuming that each particle lies along the hyperbola with 
< a/T >= 7 then each takes up a rapidity length 6 = Consequently all the n particles 

take up the rapidity length Aj/ = n6 = nm±/'j and the total area is = n^m\/^y. 

Since the particles are produced around the average hyperbola, we expect that this result 
may be modified by a constant, C 2 , of order unity giving 


hA ~ 


bc2m\rA 


Ay 


Ay 


( 20 ) 


where we have introduced C 3 = bm\c 2 - In this way we obtain from Eq.([I^ and Eq.(^0D a 
description of the grand partition function in terms of the multiplicity n. In the approxima¬ 
tion that n is large, i.e. for large rapidity intervals Ay, we can write the partition function 
in terms of two parameters ci and C 3 as 


Z ~ 


{ciAyY 


n\ 


exp 



( 21 ) 


We will comment further on the parameters Ci and C 3 when we investigate to what extent 
the partition function in Eq.(|2T|) describes the particle production in the Lund model. 


3.3 The partition function in the Gaussian approximation 

We now investigate the grand partition function in the limit where the number of particles 
is large, but the density is low (as in an ordinary gas). In this case we expect that the 
grand partition function can be approximated by the maximal term in the sum. To find the 
multiplicity for which the partition function is maximal we first define by writing Eq.(|2T|) 
as 

= exp . ( 22 ) 

If we treat n as a continuous variable we can expand <I>„ in a Taylor series as 

<h(n) ~ <I)(n) -|- (n — n) <h'(n) -|- ^"{n) . (23) 

Choosing n such that $'(n) = 0, we evidently have a Gaussian approximation for Zn 

Zn ~ exp $(rr) exp (24) 

where the variance, V, is given by U = l/<h"(n). It is straightforward to obtain expressions 
for both n and V if we use Stirlings approximation for the factorial in 4>(n). We find 

n = 


8 


^ In(^) 

2 c 3 n 


—\ -1 




V 


(25) 










Notice that, since C 3 is positive, this implies that the variance of the distribution is less than 
the mean and the distribution is therefore narrower than a Poissonian. If we now introduce 
the density of particles in the rapidity volume, R = n/Ay, then 

<I)(n) = {R + C 2 ,R^)Ay . (26) 

For large n we can approximate the grand partition function a.s Z ^ Zn and so 



with 

ciji = R c^R? ■ (28) 

The grand canonical partition function, for a gas is related to the pressure, P, temperature, 
T and volume, ln(s/so), of the gas via 


n 

P 


-kT InZ 

dQ 

dln(s/so) 


(29) 


where k is Boltzmann’s constant. For the partition function in Eg .(1271) we obtain the fol¬ 
lowing equation of state for the gas 


P = kT{R + csP^) . (30) 

Our expansion thus corresponds to the hrst two terms in the virial expansion in the particle 
density of the gas. We note that the equation of state in Eq.(^) is similar to that of a 
Van der Waal’s gas. For particles with zero (transverse) mass we have C 3 = 0. In this case 
a particle does not take up any volume in rapidity and Eq.(^) reduces to the equation of 
state for an ideal gas. 


4 The vertex distributions 

The partition function is related to the multiplicity distribution, P„, since 

Pn=A , (31) 

In the remaining sections we shall use this relationship to further study our partition function. 
We begin here with a study of the vertices produced in the string fragmentation. These 
vertices are strongly ordered in rapidity and thus satisfy one of the assumptions used to 
derive our partition function. This is only an approximation in the case of the particles. Of 
course, the number of vertices corresponds directly to the number of primary particles. 

In what follows we outline a simple model in which all particles have the same mass 
(m = 0.8 GeV) and there is no transverse momenta. The effects of relaxing those constraints 
will be considered in the next section where we return to the particles. 
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Figure 5: The logarithm of the distribution of the rapidity distance P{5y) between adjacent vertices 
as predicted by the Lund model (Eq.(|^)) for a fixed mass, m = 0.8 GeV, but different values of 
the Lund model parameters a and b. The upper plot shows fixed a = 0.3 and b = 0.4, 0.58, 1.6. 
The lower plot shows fixed b = 0.58 and a = 0.1, 0.3, 0.5. 


4.1 The distribution in rapidity 


We begin by studying the separation between neighbouring vertices. In the Lund model the 
distribution of such separations for a hxed mass, m, is given by 


P{6y) = N 


dT 



z 



/r + 

[r{i-z))) 


(32) 


The logarithm of this distribution is plotted in Fig.(^) for various values of the Lund model 
parameters a and b. We see from this hgure that there are two main characteristics of the 
distribution. The hrst is an effective minimum separation between vertices which increases 
with but is independent of a. Physically this separation arises because two vertices 
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cannot be very close together in rapidity if they mnst prodnce a massive particle. The 
second characteristic is an exponential fall off for large separations, which depends only 
on the parameter a. 

We can consider a simple model which reprodnces the above featnres very well. In this 
model the rapidity region is divided np into a series of N eqnal bins of size 5?/bin- The effective 
minimnm separation between vertices can now be taken into acconnt by demanding that no 
bin may contain more than a single vertex. Each bin is assigned a probability p to contain 
a vertex and a probability 1 — p to be empty. This allows ns to compnte the probability of 
a separation, (5p, between two vertices. If 5y is discretised as 5y = mJpbin with n an integer 
then the probability of snch a separation is given by a geometric series 


P{5y) 


p{1-pY ^ (n = 1,2 ,...) 

—^exp(-,3^!,) 

(1 -P) 


(33) 


with [3 = — ln(l — p)/5?/bin- We see that large 5y separations are exponentially snppressed. 
The two main featnres of Fig.(|^) are thns very well reprodnced by this simple model, which 
corresponds to distribnting the vertices according to a binomial distribntion (appendix ^). 

We can investigate the accnracy of the binomial approximation nsing the JETSET Monte 
Carlo (for consistency we nse a hxed mass (m = 0.8 GeV) and have no transverse momentnm 
generation). Here we generate 2-jet (qq) events and analyse the distribntion of vertices within 
a rapidity range. Ay. The energy is chosen to be snfhciently large in order to avoid edge 
effects from the q and q fragmentation contaminating the Ay region. The mean, (n), and 
the variance, V, of the resnlting mnltiplicity distribntions are nsed to calcnlate the binomial 
parameters N and p, as detailed in appendix We will see later that binomial distribntions 
with these N and p valnes do indeed reprodnce the mnltiplicity distribntions very well. Fignre 
(|) shows the resnlts as a fnnction of Ay, for varions valnes of the parameter b. We see that 
for large rapidity volnmes {Ay > 5 nnits) the binomial assnmptions seem to work very well. 
That is to say, the observed p parameter is effectively constant as a fnnction of Ay, whilst the 
parameter N is linear with Ay. This corresponds to a constant bin size (5pbin- (As expected 
the bin size is fonnd to be proportional to bm?.) 

The behavionr of the effective JETSET N and p parameters at small valnes of Ay can 
easily be nnderstood. When Ay becomes smaller than the normal bin size, Jpbin, we have 
only one bin which is now of size Ay. In Fig.(|]) all of the N cnrves indeed tend to the limit 
N = 1. Meanwhile the observed p valne is the probability to hnd a vertex in this single bin 


Pobs = P 


^Pbin 


{Ay < Jpbin) 


(34) 


Thns the observed p becomes linear with Ay when Ay is smaller than the bin size. This 
effect can be seen in Fig.(^. Between the limits of large and small Ay the behavionr of N 
and p are not so well determined. Here correlations between closely spaced vertices will play 
a role. 

We are now in a position to tnrn onr attention back to the partition fnnction. This formnla 
shonld also generate a good description of the mnltiplicity distribntion for the vertices. We 
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Ay 


Figure 6: The values of N (upper plot) and p (lower plot) for the vertex distributions produced by 
JETSET (for a fixed mass and no transverse momentum generation) as a function of the rapidity 
volume, Ay. We show the results for fixed a = 0.3 and b = 0.4, 0.6, 0.8 and 1. 


have 


Pn = Co 


icAyy 


n\ 


exp 


—bm? C2U^ 

Ay 


(35) 


Here cq is a normalisation parameter and so is determined in terms of the remaining param¬ 
eters. We can relate the parameters ci and C 2 to the parameters N and p of the binomial 
distribution. The procedure is explained in detail in appendix For large Ay, we obtain 


Cl 


C2 


Np 

Ay 


exp 

Ay 


p 


.(1 - t ). 


2bm‘^N{l — p) 


(36) 


In Fig.(0) we show the values of ci and C 2 which we obtain from our JETSET multiplicity 
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distributions. For large rapidity volumes, A?/, they tend to constant values. We noted in 
section 1^ that it is also possible to approximate Eq. (|35[) using a Gaussian distribution (with 
the appropriate mean and variance). If we express the mean and variance of Eg. (p5D in terms 
of N and p and solve for ci and C 2 , then we obtain the same expressions as Eg. (|36D . We 
note, however, that in the case of a Gaussian distribution one has a symmetric distribution. 
This is not true of either Eq.(^^ or the binomial distribution since they both contain a term 
n\ in the denominator. 

Finally in Fig.(|D we demonstrate how well the binomial and Eg. (|35D reproduce the ob¬ 
served multiplicity distribution. We show three curves hrstly the JETSET multiplicity dis¬ 
tribution, secondly that obtained from the binomial distribution and hnally the distribution 
obtained from Eg. ([5^) . At Ay = 5 we see very good agreement and it is difficult to dis¬ 
tinguish the different curves whilst at Ay = 10 all of the curves lie on top of each other. 
We thus see that the vertex multiplicity distributions produced by JETSET do indeed agree 
very well with our simple expression for the partition function, Zn- 


4.2 Distribution in proper time 


So far we have discussed the distribution of the vertices in terms of the rapidity, y. If px 
is neglected then the position of the vertices is specihed by one further variable F, which is 
related to the proper time of the vertex. In this section we will investigate how the vertices 
are distributed in F. As we discussed in section 13.11, we have for the vertices that 


F(F) oc F“exp(-6F) 


(37) 


which has a mean (F) = (1 -I- a)/h. Equation (^) is, however, an inclusive distribution. If 
we examine vertices within a rapidity range. Ay < 2, then we hnd that they are correlated. 
This means, for example, that if a vertex has a large F value then nearby vertices are also 
likely to have large F values. 

We now examine how the vertices are distributed in F inside a rapidity range Ay for 
various multiplicities, n. Motivated by the hnite energy vertex distribution, Ff(F), which we 
considered earlier in Eq. ( 0 ). we parameterize the distributions as 

pagfj(n,Aj/) 

( 38 ) 

Note that taking the weighted average of Hn should reproduce the inclusive distribution of 
Eq. (^TD . In Fig.(^ we show distributions in F obtained from JETSET for Ay = 6. Each 
plot in the hgure is for a different number of vertices, n, together with the corresponding 
£t according to Eq. (|38|) . Here the parameters Oefr and Seg have been fitted for each different 
n value. We see that one can find values of aeff(n. Ay) and Sefr(n, Ay) for which a very 
reasonable description of the F distributions is obtained. We note that the large F behaviour 
is determined only by the Lund parameter b and not by n or Ay. Thus it is only dependent 
on the scale for the area law suppression. Next we examine the dependence of both 
and Seff on the multiplicity and the rapidity interval. We have carried out hts to the F 
distribution obtained from JETSET for a set of values of n and Ay. We hnd that both of 
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Figure 7: The values of ci (upper plot) and C 2 (lower plot) for the vertex distributions produced by 
JETSET (for a hxed mass and no transverse momentum generation) as a function of the rapidity 
volume, Ay. We show the results for fixed a = 0.3 and b = 0.4, 0.6, 0.8 and 1. 

these functions depend only on the ratio R = n/Ay. This can be seen clearly in Fig.(^) 
where we plot the results of our hts for three different values of Ay, as a function of the 
density R. 

This completes our study of the distribution of vertices produced in the Lund model 
of fragmentation. We can summarize our hndings as follows. In rapidity the vertices are 
approximately distributed according to the partition function, whilst in proper time they 
are distributed according to Eq. (|38D . Importantly we hnd that the large F behaviour of the 
distribution in F is determined only by the area law. We also hnd that the functions aes and 
Seff only depend on the density of vertices, R, which is itself the important quantity in the 
equation of state for the gas in rapidity. 
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Figure 8: The multiplicity distributions for the vertices for various values of the rapidity volume 
Ay as obtained from JETSET (solid curve), the binomial distribution (medium grey curve) and 
Eq.d^) (light grey curve). 

5 The particle distributions 

For primary particles the mean multiplicity corresponds to the mean number of vertices. The 
effects of going over from vertices to particles essentially means some smearing in rapidity. 
Thus the rapidity ordering assumed for Eq. (|^) will no longer be true. However, for large 
Ay it should still be a good approximation. The rapidity of a particle is distributed around 
the average rapidity of the two vertices from which the particle stems with a width of about 
one unit of rapidity. Therefore the particle multiplicity distribution for a hnite rapidity 
interval Ay will have the same average as the vertices but a larger width. To understand 
this effect we return to our simple binomial model. As before we divide the rapidity range 
into N equal bins with the probability p to contain a vertex. Now we further assume that 
the presence of a vertex in any bin results in a particle in one of the two neighbouring bins 
with probability q or in the original bin with the remaining probability 1 — 2q. In order to see 
how this smearing affects the mean and the variance we compute the generating function. 
The generating function for the original binomial distribution is given by 

G{z) = [l+p{z-l)f . (39) 

Some straightforward algebra then shows that the generating function for the above particle 
distribution is given by 

G{z) = [1 +p{z - 1)](^“^)[1 +p{z - 1) +p‘^{z - l^qil - q)f ■ (40) 

Thus two factors of (1 +p{z — 1)) have been modihed. The mean is unchanged and equal to 
Np, but the variance is increased from V = {n){l — p) to 

V = {n)[l —p +Apq{l —q)/N] . (41) 

This distribution can be rather well approximated by another binomial distribution with the 
same mean and variance. This corresponds to effective p and N values 

Peff = t[1 - 4g(l - q)lN] 
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Figure 9: The distributions Ay) obtained from JETSET for the default values of the Lund 

model parameters (a = 0.3 and b = 0.58). Also shown are the continuous curves obtained from our 
fits based on Eq.(^). In this example Ay = 6 and n = 3... 10 . The corresponding R values are 
shown on each plot. 




Np 

PeS 


(42) 


Thus we see how a larger spread of the particles around the vertices (a larger g-value) 
corresponds to a larger width and a smaller effective p-value. From Eq. (|4^ pefr must be 
larger than zero, but if we had allowed for a spread beyond the nearest bin then negative 
values of Peff would be possible. This corresponds to a negative binomial distribution. Since 
the bin width is of the order of bm^ the spread is certainly beyond neighbouring bins in the 
case of pion production. 

We have investigated various cases of final state production. The multiplicity distributions 
can still be well approximated by binomial distributions with constant p-values for large 
rapidity intervals. In Fig. (|lT|) we show N and p as a function of Ay for the various cases. 
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Figure 10: The functions aeff(ii) (upper plot) and Sefr(i?) (lower plot) versus R. We show the 
results for three different rapidity ranges Ay = 4 (triangles), Ay = 5 (circles), and Ay = 6 (stars). 


For a situation with only a single stable hadron, assumed to have the mass m = 0.8 GeV, 
and no transverse momentum generation, the result is as expected. Comparing the multi¬ 
plicity distribution with the distribution for the vertices, we hnd that p is decreased and N 
is increased. The product of N and p is however the same for the two distributions. 

If we include the standard mixture of different hadron masses p is further reduced. We 
obtain in this case a distribution that is very close to a Poissonian. Thus, as expected, the 
width of the multiplicity distribution greatly increases when light pions are produced. 

Including transverse momentum generation increases p to positive values as shown in the 
hgure. The transverse mass of the pions is thus, in the case of the standard mixture of 
hadrons, not small enough to give a negative p-value. 

Finally, if we include the decays of unstable particles and analyse the hnal charged particles 
then the width increases substantially and p becomes negative (corresponding to a negative 
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Figure 11: The values of N (upper plot) and p (lower plot) for the particle multiplicity distributions 
produced by JETSET. The four cases shown correspond to: (a) a single mass, no decays or p±-, 
(b) complete mass spectra, but no decays or p±; (c) complete mass spectra, but no decays; and (d) 
complete mass spectra, charged hnal state particles. 


binomial distribution). Including the final uncharged particles in the analysis results in an 
even more negative p-value. 

We can summarize our hndings as follows. The width of the multiplicity distribution is 
very sensitive to the mass spectrum of the produced particles. Using default JETSET the 
average transverse mass is large enough to give a binomial multiplicity distribution. In this 
case the negative binomial distribution for the hnal state stems from the increased width 
due to decays. 
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6 Conclusions 


Inspired by the Feynman-Wilson gas analogy we have derived an explicit form for the grand 
partition function of the Lund fragmentation model. This partition function is described 
in terms of the multiplicity n. In particular, we derive an equation of state for the gas, 
corresponding to the first two terms in the virial expansion in the particle density. 

The partition function is derived in the approximation that the particles are ordered in 
rapidity. This is true for the string break-up vertices and the number of vertices corresponds 
to the number of particles. Therefore, we have investigated the properties of the partition 
function using the vertices. For large rapidity intervals, we find that the average and the 
fluctuations of the multiplicity of vertices are described by the partition function. 

The partition function gives a multiplicity distribution which is close to a binomial dis¬ 
tribution. We hnd that the average transverse mass of the produced particles is sufficiently 
large to get a reasonable description from the approximation that the particles are ordered 
in rapidity. Thus the multiplicity distribution of the particles stemming from the string is 
described by an ordinary binomial. It is the decays of the unstable particles that results in 
a negative binomial distribution for the number of final charged particles. 

The distribution of the vertices for different rapidity volumes and different multiplicities 
has also been investigated in terms of the proper-time. We find that the behaviour for large 
proper-times is determined only by the area-law and is independent of both the volume 
and the multiplicity. For smaller proper-times the distribution is described by a simple 
parametrisation. We hnd that the important quantity for the parametrisation is the density 
of vertices in rapidity, which in turn is described by the equation of state for the gas. 
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A The binomial and negative binomial distributions 

The binomial distribution is dehned by 

P(n)= ( ^ ) (43) 

The average, (n), and the variance V = {v?) — (n)^ of this distribution are related to N and 
p via 

(n) = Np 

V = Np(l - p) . (44) 
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The binomial distributions form a family of distributions depending on the values of N 
and p. In the limit p —0 for constant (n), the distribution becomes a Poisson distribution. 
It is also possible to continue the expressions in Eq.(^) to negative p-values, which for 
constant (n) = Np implies also a negative N. In this case the distribution becomes a 
negative binomial distribution. Such a distribution is conventionally written in the form 


Pk{n) = ^ + ^ ^ ^ ^ p\l-p) 

where 

P = (pO)) 

k = -N {N <0) . 


(45) 


(46) 


Note that the relationships of Eq. (|4^) for the average and the variance remain true, even 
when p and N are both negative. Negative p-values correspond to distributions which are 
wider than a Poissonian. Thus the negative binomial distributions belong to the same larger 
family as the (ordinary) binomial and Poisson distributions. Within this family the width 
can vary from zero to infinity. Ordinary and negative binomials correspond to V smaller and 
larger than (n) respectively, with the Poisson distribution as the limiting case in between. 


B The binomial approximation of the partition func¬ 
tion 


Our aim is to determine the ci and C 2 parameters of Eq. (|35[) from the N and p parameters 
of a binomial distribution. We begin by using Stirlings approximation to write the binomial 
distribution as 


ln(F) = 


InN 


+ N \nN -{N-n) \n{N - n) 


n 


ln(iV — n) 


+ iVln(l — p) + n In 


p 


1 — p 


— ln(?7,!) 


(47) 


whilst the distribution in Eq.(p5D can be written as 

ln(F„) = In(co) - + n ln(ciA 2 /) - ln(n!). (48) 

Ay 

We now express n as n = (n) + x, where {n) is the mean. Next we subtract Eq.(^7D from 
Eq.(^8l) and expand around a; = 0 up to terms of order Equating the series coefficients 
to zero determines the parameters ci and C 2 in terms of N and p. We obtain 


CiAy 

bm^C2 

Ay 


Np exp 


1 — 2p + 2Np — 2Np^ 


2N{l-py 


2N{l-p) - 1 
4A2(1 -p)2 


(49) 
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Which for large A?/ can be simplified to 


Cl A?/ = Np exp 
brn^C2 1 


P 


_(1 - p )_ 


Ay 2N{1 — p) 

If we insert these expressions into Eg. (|5BD then we finally obtain 




{NpY 


n\ 


■ exp 


+ n2Np 
2N{l-p) 


(50) 


(51) 
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